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We present an experimental implementation of the coined discrete time quantum walk on a square 
using a three qubit liquid state nuclear magnetic resonance (NMR) quantum information processor 
(QIP). Contrary to its classical counterpart, we observe complete interference after certain steps 
and a periodicity in the evolution. Complete state tomography has been performed for each of 
the eight steps making a full period. The results have extremely high fidelity with the expected 
states and show clearly the effects of quantum interference in the walk. We also show and discuss 
the importance of choosing a molecule with a natural Hamiltonian well suited to NMR QIP by 
implementing the same algorithm on a second molecule. Finally, we show experimentally that 
decoherence after each step makes the statistics of the quantum walk tend to that of the classical 
random walk. 
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I. INTRODUCTION 

The idea of exploiting the quantum mechanical be- 
haviour of a device to gain power in simulating quan- 
tum systems was first introduced by Richard Feynman 
0. The field of quantum computing has since grown 
enormously with the discovery of two algorithmic pillars: 
Shor's factoring algorithm 2] and Grover's search algo- 
rithm Both of these demonstrate a clear speed-up 
over their classical counterparts. Following in this path, 
many other quantum algorithms have been developed 
that provide a speed-up 0, 0, IB| ■ A more recent addi- 
tion to the family of quantum algorithms which demon- 
strate an exponential speed-up, are those based on the 
quantum random walk - the quantum version of the suc- 
cessful classical random walk 7] . 

There is however, a need to explore more than the 
simple computational properties of the algorithms. They 
must also be experimentally tested in real devices and 
their relative ease of implementation compared and con- 
sidered. In particular, in QIP devices where we are 
controlling the natural Hamiltonian, it is important to 
choose a system where the Hamiltonian is amenable 
to automatic and systematic control. This can be ex- 
plored by implementing the same algorithm in different 
molecules and contrasting the performance. Although 
many different implementation schemes have been pro- 
posed for the quantum random walk algorithm, using 
for example trapped ions 8J, an optical lattice 9], cav- 
ity QED 10], or an optical cavity [ll|, these have not 
been tested. The only experimental test of a quantum 
walk is the continuous time version of a quantum walk 
on a square using a two qubit nuclear magnetic resonance 
(NMR) quantum information processor (QIP) [13. This 
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work showed the contrast between a classical and quan- 
tum random walk and showed the influence of entan- 
glement on the probability distribution of the quantum 
walk. Here, we present an experimental proof of principle 
experiment of a discrete time quantum walk on a square. 
The effects of decoherence on the quantum random walk 
has been investigated by several authors and indeed, it 
may offer some benefits ■ Therefore, we also ex- 

plored the quantum to classical transition of our walk 
under the addition of decoherence to the coin. Further- 
more, we compared and contrasted two different control 
schemes and molecules by implementing the algorithm 
on two molecules. 



II. QUANTUM RANDOM WALKS 

In the development of deterministic classical random- 
ized algorithms, the methods of Markov chains and ran- 
dom walks have played a fundamental role These 
algorithms can be divided into two categories: contin- 
uous time random walks when the walker has a proba- 
bility per unit time to make a move; and, discrete time 
random walks where the walker moves at defined time- 
steps. Since these processes are stochastic, it is not sur- 
prising that they have quantum counterparts. The quan- 
tum versions however, show remarkable differences with 
their classical analogues. The continuous time quantum 
walk (CTQW)0| has been shown to provide an expo- 
nential speed-up in propagation through a grap h [3, 113 • 
The discrete time quantum walk (DTQW) [l^ plays an 
important role in the speed up of a quantum algorithm 
design for spatial searching [I^ |20, |2j| . 

One step of a classical discrete time random walk on 
a circle with n nodes, denoted by {0, . . . ,n — 1}, is per- 
formed by repetition of the following two steps: (1) the 
walker first fiips a coin and then (2) moves either clock- 
wise or counterclockwise depending on the outcome of 
the coin toss. 
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FIG. 1: Comparison of the dynamics of the classical (left) and quantum (right) random walk on a square for three steps. H 
or T represents the state of the coin and the number, the position of the walker at one of the four nodes of the square, p is the 
probability of each classical state and A is the probability amplitude for the quantum state. Part (a) for each step is the coin 
flip and part (b) the movement around the square. In both cases the walker starts at node with the coin in the heads state. 
After one step he has a fifty percent probability of being at either of nodes 1 or 3. Then, in the second step he goes to either 
or 2 with fifty percent probability. In the third step however, the two types of walk diverge. The classical walk continues to 
oscillate and the probability remains spread out. In the quantum walk on the other hand, the probabilities interfere and cancel 
out leaving all the probability in one corner after three steps. 



If we perform a quantum mechanical treatment of the 
situation, we can label the nodes with a mutually or- 
thonormal set of state vectors {|i)}"Jo^. The coined 
DTQW on the circle can be seen as "quantumly" flip- 
ping a coin degree of freedom using a unitary operation 
and then coherently moving the walker position degree 
of freedom clockwise, or counterclockwise, conditioned 
on the state of the coin 22]. For a Hadamard walk, the 
coin flipping operation is simply the Hadamard gate de- 
scribed by the matrix, 



Now, the conditional shift operator is defined as 



S\Hm = \H)\iQl) 
S\Tm = |T)|*©1), 



(1) 



(2) 
(3) 



where © and are understood to be addition and sub- 
traction modulo n and \H) and |T) describe the two basis 
states of the coin. Therefore, if the walker is in position 
he will move clockwise to the position |i G 1) if the 
coin in the state or counterclockwise to \i © 1) if the 
coin in the state \T). We can write this operator as 

S = Y,{\H){H\®\^Ql){^\ + \T){T\®\^®l){^\).{A) 

i=0 

Then, one step of the DTQW is defined as applying 
the operator 



On a circle, this type of algorithm shows destructive 
interference effects and a probability distribution that is 
periodic in time. The contrasting dynamics for the clas- 
sical and quantum random walks are shown in Fig. ^ 
As opposed to the classical walk where the probability is 
always spread out, the quantum walk has steps where the 
probability amplitudes interfere such that all the prob- 
ability comes back to one node. Furthermore, this walk 
is periodic in that after eight steps, the corresponding 
propagator is equal to the identity and the system comes 
back to its original state. 
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W = S{H®t). 



(5) 



FIG. 2: Logical labeling of the nodes on which we imple- 
mented the DTQW. With this labeling, flipping the flrst qubit 
corresponds to a horizontal move and flipping the second 
qubit, a vertical move. 

In our experimental setup we have three qubits avail- 
able, which allows one qubit to describe the coin state 
and two for the position state. Thus, we have n = 4, and 
we are performing a discrete quantum walk on a square. 
The shift operator defined in Eq. 0]would require a com- 
plicated quantum circuit involving a Toffoli gate. We 
can simplify the circuit required by using a shifting op- 
erator that moves the walker along a direction vector. 
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i.e. horizontally or vertically (this also is analogous to 
the random walk on the hypercube ^\). Therefore, if 
we label the corners of the square as shown in Fig El the 
shift operator on the three qubit register becomes, 

S ^ PhX'^ + PtX^ 

= {PhX^ + Pt){PtX^ + Ph) 

= {X^Cnof-^X^)Cnof-^, (6) 

with X denoting the standard ax Pauli matrix, Ph/t S'^e 
the projectors on the two coin states, and the superscript 
indicates on which of the qubits the action is performed. 
Here, it is understood that the first qubit represents the 
coin and the second and third, the position register. The 
resulting probabilities for each step are shown in Table 
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TABLE L Probability to be in each of the corner states as 
denoted in Fig. |5| While in the classical random walk the 
probability always remains spread out between two corners, 
in the quantum random walk all the probability returns to 
one corner at certain time steps. 



III. LIQUID STATE NMR QUANTUM 
INFORMATION PROCESSING 

A. The basic principles 

A liquid state NMR QIP consists of an ensemble of 
roughly 10^° identical molecules dissolved in a liquid sol- 
vent. Due to the fast tumbling motion of the molecules, 
they are essentially decoupled from each other; ideally 
all the molecules have the same evolution. We can think 
of the quantum register made of qubits that correspond 
to the spin ^ nuclei within each molecule. The sample 
is placed in a strong homogeneous magnetic field which 
provides the quantization axis and causes the spins to 
precess around the axis of the field. It is possible to 
implement single qubit gates using radio-frequency (r.f.) 
pulses resonant with the precession frequency, which can 
effect a rotation about any axis orthogonal to the axis of 
the external field. Two qubit gates are effected through 



the coupling from the natural Hamiltonian, which pro- 
duces a controlled phase gatej24|. 

If the molecule used contains n distinguishable nuclei 
and the magnetic field is aligned along the z-axis, then 
the system Hamiltonian is approximated by, 

n 

n = T:^viZ, + '!^^j,jZ^® Zj, (7) 

i—l i>j 

where Vi is the Larmor frequency of spin i in Hz, Jy is the 
coupling strength between spin i and j in Hz and Z in the 
conventional Pauli operator ct^. The interaction part of 
the Hamiltonian can be approximated to the above Ising 
form (weak coupling regime or secular approximation) 
only if the difference between any two nuclei Larmor fre- 
quencies is much greater than the coupling between the 
nuclei. We can also turn off the coupling between any 
two spins as needed by applying refocussing r.f. pulses. 

B. Implementing dephasing in NMR 

We can apply a controllable amount of decoherence to 
selected spins using gradient techniques in NMR. Con- 
sider only one nucleus with state p and suppose we work 
in the rotating frame of that spin. On a NMR spectrome- 
ter, it is possible to apply a gradient to the external mag- 
netic field. During the time that the gradient is applied, 
the spins will precess at different frequencies depending 
on their position in the sample. The state of the ensem- 
ble will then be given by an average over the observable 
sample, 

where 2a is the length of the sample, t is the interval 
of time the gradient is being applied and a' ~ a/h and 
7 = V jBzi the gyro-magnetic ratio of the nucleus. If we 
compute the integral, it can be shown that, 

p' = (l-p)p + pZpZ, 

p=i('l-^i-sin(a'7to)), (9) 

which is the exact form of a z-dephasing decoherence. 
The amount of dephasing can be controlled by the 
strength and time of the gradient pulse. Particular spins 
can be protected from the applied decoherence by apply- 
ing a 180 degree rotation and applying a second gradient 
of the same strength and time. This second gradient will 
reverse the dephasing of the rotated spins and double it 
on the spins that were not rotated. 

IV. THE EXPERIMENT 

We implemented the quantum walk algorithm on 
two molecules: trans-crotonic acid and trichloroethylne 
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(TCE). This allowed us to compare the quality of two 
different methods of control and the merits of the two 
molecules. 



A. Implementation on crotonic acid 

The seven qubit molecule trans-crotonic acid (four car- 
bons, two hydrogens and one methyl group) has been 
used in experimental demonstrations of quantum algo- 
rithms, such as quantum error correction |2^, and 
quantum simulations |27j. In this experiment, we used 
the carbon back-bone of labeled trans-crotonic acid in 
a solution of deuterated acetone. The hydrogen nuclei 
were decoupled using standard heteronuclear decoupling 
techniques 28] . We used C3 as the coin and C2 and C4 
as the position register (see Fig. OJ. Ci was used as 
a labeling spin to ease the creation of the initial state. 
On a Bruker DRX Avance 600 NMR spectrometer, the 
molecule has the Hamiltonian parameters shown in Fig. 
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FIG. 3: Molecular structure of trans-crotonic acid and its 
Hamiltonian parameters. The chemical shifts are given as 
the diagonal elements and the coupling strength (Hz) by the 
off-diagonal elements. Note that since the darkly shaded un- 
labeled nuclei are oxygen whose natural abundance of ^^O 
with spin is close to 100 %. Therefore, the two oxygen nu- 
clei do not couple with the rest of the molecule and can be 
ignored. Lightly shaded unlabeled nuclei are hydrogen which 
were decoupled during the experiment. 



Since our system is homonuclear, the control of indi- 
vidual qubits is achieved through soft gaussian-likc r.f. 
pulses at the Larmor frequency of the target nucleus. 
The length of the soft pulses is of the order of the in- 
verse of the smallest chemical shift difference with the 
other nuclei. In our experiment the length of the selec- 
tive pulses on Ci and €2,03,04 were 192^s and 704/is 
respectively. 



1. Initial state preparation 

The experiment required the initial state 

p,„^*(E) |000)(000| = *(]1 +Z){t + Z){t + Z). (10) 

We created the labeled pseudo-pure state ZOOO (using 
the notation C1C2C3C4) following the spatial averaging 
technique elaborated in p^. 

2. Pulse sequence implementation 

The unitary of one step of the DTQW from Eq^jwas 
translated to a sequence of pulses and coupling gates as 
shown in Fig. ^ Although many pulse sequences are 
possible through the use of commutation rules, this par- 
ticular one was designed to be the most efficient due to 
the cancellations possible during multiple step sequences. 
Moreover, the ZZ gates are achieved simultaneously, 
which shortens the overall pulse sequence, thus reduc- 
ing decoherence effects. Commutation rules were also 
used to cancel pulses between the final step and the read- 
out pulses. The ideal pulse sequence of rotations and 
couplings was then input into a pulse sequence compiler 
which numerically optimized the timing and phases of the 
pulses js^]. The compiler pre-simulates the selective r.f. 
pulses using an efficient pair- wise simulation and then de- 
composes the simulated unitary into phase and coupling 
errors sandwiching the ideal selective pulse. These errors 
can then be taken into account by the refocussing scheme 
and phase of the pulses, so that the overall unitary is as 
close to the desired one as possible. 




FIG. 4: NMR pulse sequence representing one step of DTQW. 
The notation i?f means a rotation of an angle 9 around the 
axis i. Refocussing pulses are not shown. Since each nucleus 
is tracked in its own rotating frame, rotations about the z- 
axis are implemented instantaneously through a change of 
reference frame. 

Since we are concerned with the final state of only three 
qubits in this experiment, complete state tomography is 
still feasible. On a three qubit system in NMR, only seven 
different readout pulses are required to rotate each term 
of the density matrix into observable simple single coher- 
ences [3^ . And, since we were operating on a homonu- 
clear system, observing the signal from all spins in one 
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experiment was possible, with some post-processing to 
adjust for the correct phase of each individual rotating 
frame. The coupling between the labeling spin Ci and 
the other three qubits is resolvable and so the presence of 
the labeling spin does not interfere with the tomography 
of C2, C3 and C4. 

3. Experimental results 

For the state tomography each of the peaks in the spec- 
tra were fitted using absorption and dispersion Lorenzian 
peaks. The full density matrix was then reconstructed 
by appropriately summing up the corresponding Pauli 
terms. Where two experiments gave values for the same 
density matrix terms, the values were simply averaged. 
As we observed only C2 , C3 and C4 , the term ZIII could 
not be determined. A suitable amount of that term was 
subsequently added to the density matrix so as to make 
the initial state as close to ZOOO as possible. This amount 
was then kept constant for the density matrix reconstruc- 
tion in subsequent experiments. 

To quantify the success of our experiments, we com- 
puted the fidelity of the experimental density matrix to 
both the ideal and simulated results. In NMR, all states 
are nearly completely mixed and the fidelity measure in- 
troduced in 1^ is appropriate. We can compare one 
density matrix to another using the following formula. 



We made two comparisons. Firstly we compare the ex- 
perimentally determined density matrix to the theoreti- 
cally expected result. The theoretical result is achieved 
by multiplying the ideal initial state by the ideal propa- 
gator. To investigate how well we understand our control 
of the system, we also compare the fidelity of the results 
from a simulation of the experiment to the theoretical 
result. 

The fidelities of simulated and experimental results are 
compared in Table ^ The loss of fidelity in our experi- 
ment, over and above that of the simulated control errors 
is explained from three sources which are not taken into 
account by the simulation. We have losses from T2 re- 
laxation. Although our pulse sequence is short compared 
with the T2 relaxation times, during the quantum walk 
algorithm, the state is often in high coherences, which 
decay much faster than the simple T2 time. Inhomo- 
geneities in the strong magnetic field also cause extra re- 
laxation and dephasing. Further losses come from inho- 
mogeneities of the r.f. field used to implement rotations 
and pulse angle mis-calibration. 

4. Addition of decoherence on the coin 

In a subsequent experiment, we added dephasing de- 
coherence to the entire qubit register using the technique 
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TABLE IL Fidelities (in percent) of experimental and simu- 
lated results. The first column gives the fidelity of experimen- 
tal density matrix determined from the tomography, with re- 
spect to the theoretical expected density matrix. The second 
column gives the fidelity of the simulation results. Errors are 
estimated from the fitting procedure. Note that since com- 
puter simulation of the spatial averaging that occurs during 
the pseudo-pure preparation is difficult and inaccurate, the 
initial state for the simulation was the experimental pseudo- 
pure state determined from the tomography. The fluctuations 
observed in the fidelity come from uncertainties in the fit and 
instabilities in the spectrometer over the course of the exper- 
iment. 

described in section IIII Bl We expect the behaviour of 
the quantum walk should converge to the classical walk 
as the decoherence becomes complete after each step. To 
demonstrate this claim experimentally, we implemented 
the quantum random walk for four steps, adding decoher- 
ence of a certain strength between each step of the walk. 
The diff'erences between quantum and classical walk is 
manifested in the different probabilities of being in each 
of the corners after each step. The results are shown 
in table IIIII for gradient strengths corresponding to no, 
partial and full decoherence. 

The divergence between the classical and quantum 
walk shows most clearly in steps three and four. Whether 
the walk is classical or quantum, steps one and two 
yield the same measurement probabilities for the posi- 
tion (however the quantum version with decoherence will 
have coherent superposition states). Analyzing the data 
from step 3 and 4, one can see that the quantum inter- 
ference present so clearly in the quantum walk with no 
decoherence, is less obvious as the amount of decoher- 
ence increases. Instead of the probability all collecting in 
one corner it remains spread out between two opposite 
corners - the same as in the classical walk. 

The probabilities even with zero gradient strength do 
not correspond perfectly to the ideal quantum walk. We 
believe these errors come from two sources. Because the 
gradient does not commute with any pulses we were not 
able to use commutation rules to reduce the number of 
pulses during multiple step experiments. Furthermore, 
gradient methods are hampered by diffusion and multiple 
gradients may lead to a return of signal that was "erased" 
by a previous gradient. 
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TABLE III: Classical versus estimate of quantum probability to be in each corner of the square for one through four steps (c.f. 
Table HJ. The quantum results were obtained for gradient strengths corresponding to no, partial and full decoherence. The 
experimental results were obtained by reconstructing the density matrix using the same fitting software used before and then 
applying the position measurement projectors to the reconstructed density matrix. 



B. Comparison with the TCE molecule 

For comparison purposes and to show the importance 
of choosing a molecule with good characteristics in liquid 
state NMR quantum information processing, we show our 
results from our initial attempt to implement the quan- 
tum walk on the molecule tri-chloroethylene (TCE) - a 
molecule with which we have much less control due to the 
presence of strong coupling. The molecule has been used 
for some initial demonstrations of quantum algorithms 
[3^ . A diagram of the molecule and the parameters of 
its Hamiltonian are shown in Fig. [S] 
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FIG. 5: Diagram of "C labeled TCE. The chemical shifts 
and couplings are given in the table. Note that since the 
chlorine nuclei (unlabeled) have a spin of |, they have an 
electric quadrupole moment which causes them to decohere 
quickly and they have a very small coupling to the rest of the 
molecule which we can ignore in the natural Hamiltonian of 
the molecule. 



1. Pseudo-pure state preparation 

Since the TCE molecule contains only three qubits, we 
are unable to create the labeled pseudo-pure state that 
we used in the crotonic acid experiments. Instead, we 



chose to use temporal averaging and add three separate 
experiments to achieve the initial state |000). The three 
different initial states we used are, 

pi = Z(S{t + Z)®{t+Z) 

P3 = ]l(g)]l®(]l + Z) (12) 

If we add the results of these three experiments, it 
is equivalent to having performed the algorithm on the 
initial state 

Ptn = Pl+ P2+ P3 

= il + Z)^il + Z)^{t+Z) 

= |000)(000| (13) 

Since there is only one hydrogen nucleus in the 
molecule, we can use broadband hard pulses to con- 
trol it. One useful property of the TCE molecule in a 
600Mhz spectrometer, is that the J-coupling between the 
two carbons is almost exactly 10.5 times smaller than 
the difference in chemical shift (Aiy). Therefore, dur- 
ing the time for a ^ZZ coupling gate between the two 

carbons (At = ttj^ — ), the relative chemical shift evo- 
^ 2JC1C2" 

lution of C2 with respect to Ci will be 6* = 27rAz/At — 

-fAiL. ^ _l0.57r = -7r/2 mod27r. Therefore, in the refer- 

'JC-1C2 

ence frame rotating at the Larmor frequency of Ci , every 
time there is a 7r/2 coupling between the carbons, an ex- 
tra Rz is naturally performed on C2. 

The chemical shift difference between the two carbons 
is small and the coupling between them large, so selective 
pulses were impossible to achieve using the same tech- 
nique of gaussian-shaped pulses used in the crotonic acid 
experiments. These pulses would be very long (roughly 5 
ms) and the large coupling errors that would occur dur- 
ing the pulse would be difficult to refocus. Instead, it was 
possible to perform selective pulses using hard pulses and 
the chemical shift evolution. To illustrate the technique, 
we demonstrate how to perform a selective ^ rotation of 
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C2- If we use a reference frame rotating at the Larmor 
frequency of Ci, then, during a time r = the spin 
Ci will not precess while C2 will undergo a rotation of 
— 7r/2 around the z-axis. Since, is much less than 
the coupling time ^ , we can ignore the coupling be- 
tween the two carbons and refocus only the hydrogen. 
Using this selective 2-rotation combined with hard pulses 
which rotate the two carbons together, we can perform a 
^ rotation with phase (j) on only C2 as follows: 



7r/2 
0+7r/2 



11(g) i?' 



n/2 
)-7r/2 



^-^/2^V2 



0+7r/2 



11 



(14) 



Similar pulse sequences can be derived to perform a tt 
rotation on C2 and selective pulses on Ci . Because of the 
different form of selective pulses used, the pulse sequences 
were written and optimized by hand. This required a 
different pulse sequence implementation of the quantum 
walk unitary which avoided as much as possible selective 
pulses and 2;-rotations where possible. The one 2:-rotation 
used, is a natural outcome of the Ci — C2 coupling gate 
as described above. This alternative pulse sequence is 
shown in Fig. |B1 
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FIG. 6: Circuit used to implement one step of the DTQW on 
the TCE. The z-rotation on C2 occurs naturally during the 
coupling with C\. Note also that the refocusing pulses are 
not shown in this pulse sequence. 



2. Experimental Results 

Fidelity results, similar to those calculated for the cro- 
tonic acid experiments are shown in Table IIVI Clearly 
this experiment was not as successful as the implementa- 
tion on the crotonic acid molecule. There are two main 
reasons for this loss of fidelity. Firstly, the chemical shift 
difference between the two carbons is very small. Be- 
cause of this, the secular approximation no longer holds 
and thus the coupling between the two carbon spins can 
no longer be approximated by the Ising form Zc-^ ® Zc^ ■ 
Indeed, it has to take all the strong coupling terms into 
account, i.e Sci-Sc2 = Xc^^Xc^+Yc^^Yc^+Zc^^Zc^- 
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TABLE IV: Experimental and simulated fidelities (in per- 
cent) for the implementation of 8 steps of the DTQW on the 
molecule TCE. Again, simulation of the pseudo-pure state 
was not performed. 



Unfortunately, this strong coupling renders our ideal 
ZZ gates much less precise. Every coupling gate per- 
formed added XX and YY error terms which we could 
not refocus. This coupling also caused problems dur- 
ing our selective carbon rotations. Although the cou- 
pling is small, there is an unrefocusable coupling of 
^4Ai^^ ~ 4.27°. Our only way to minimize these er- 
rors was to optimize the delay times analytically and 
from numerical simulations. However, these did not cor- 
respond well to the experimentally determined optimal 
values. This point also clearly demonstrates the second 
reason for the less satisfactory results on TCE. We were 
unable to use the numerical optimization of the pulse 
sequence compiler used for crotonic acid. The compiler 
provides a systematic and reliable way to produce pulse 
sequences which implement unitaries with high fidelity 
and is clearly superior to writing and optimizing pulse se- 
quences by hand. These experiments also showed the lim- 
its of our simulator. For the crotonic acid experiments, 
where only soft pulses were used, the r.f. power applied 
changed slowly and the simulator was faithful to what r.f. 
power the spins were experiencing. In TCE, where con- 
trol was achieved only through short hard pulses, other 
effects such as phase transients enter and the spins might 
experience an r.f. field much different from the ideal 
square pulse simulated. To fully understand the issues 
surrounding hard pulse control a much more detailed 
study of the probe response must be undertaken. This 
underlines a key point: control of a more complex and 
strongly coupled system could be obtained through so- 
phisticated control techniques such as strongly modulat- 
ing pulses j3lj | ; however, it seems prudent to invest the 
effort in a wise choice of molecule. 



V. CONCLUSION 

We have presented the first experimental implementa- 
tion of a coined discrete time quantum walk. It showed a 
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clear difference with the classical coined quantum walk, 
since the DTQW possesses destructive interference and 
periodicity in its evolution. A proof of principle like this 
lays down the path to more elaborate experiments us- 
ing discrete quantum walks, such as the database search, 
walks on hypcrcube or N-nodes circle or a more profound 
study of the effect of decoherence on the walk. This paper 
also demonstrates the importance of choosing a natural 
Hamiltonian well suited to automated control in the con- 
text of quantum information processing. 
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